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STRESS ANALYSIS AND STATIC TEST OF A THREE CELL BEAM 


INTRODUCTION 


An aircraft wing consists of structural members capable of resisting 
bending, torsional and shear loads. Insofar as the stress analysis of the 
wing is concerned, the structural members resisting each of these loads 
can be considered separately. In the consideration of bending loads, wing 
structures may be typed according to the bending-load resistant material. 





Yor example, all bending material may be concentrated in the spar caps or 
distributed around the periphery of the profile. The analysis of the 
first type treats the wing covering as ineffective in bending. This is a 
reasonable conclusion, inasmuch as the sheet-metal covering will buckle at 
a very low load and the load carrying ability insofar as bending is con- 
cerned, is therefore negligible. As for the second type, the distributed 
bending material consists of spanwise stiffening members. These members 
are in the form of stringers, corrugations, and inner skins fastened to the 
covering metal. However, on the compression side of the wing the metal 
sheet covering will buckle at relatively low loads which necessitates the 
use of effective widths. Hence, the stress analysis of such a structure 
considers the skin as partially effective in bending. 

At present there is the necessity of constructing aircraft wings with 
comparatively thick metal sheet covering. This increase of covering thick- 
ness permits the wing to carry relatively high airloads before buckling 
of the skin on the compression side, thus the skin is fully effective in 
bending. The stress analysis of this type of wing is the subject of this 


report. 
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PART I 


ANALYSIS 


THREE CELL BEAM WITE SKIN FULLY EFFECIVE IN BENDING 








PART I 
ANALYSIS 
THREE CELL BEAM WITH SKIN FULLY EFFECTIVE IN BENDING 


The object of this analysis is to develop the shear flow equations 
for a three cell beam with skin fully effective in bending. The re- 
sulting equations will apply for a beam loaded with shear, bending and 
torsional loads. 

Consider a three cell beam with a general profile (e.g., any airfoil 


shape) as shown in Figure 1. 





Refering to Figure l. 


0 centroid of cross section. 


distance from centroid to shear center. 


b> 
Ni 


V = shear force in y direction (through shear center). 


Initial Assumptions: 
1. Principal axes are y and x axes, where x axis is chord line. 
2. Wing planform contains no taper. 
3. Load applied through shear center in y direction. 


The shear flow designations are as shown in Figure e. 











Figure 2 


Referring to Figure 2. 


qa, = Shear 
1, = Shear 
a, = Shear 
q. = Shear 
qa = Shear 
Lo = Shear 
a. = Value 
qe = Value 
= Valu 
qs alue 


flow at any point 
flow at any point 
flow at any point 
flow at any point 
flow at any point 
flow at any point 
of q, at point a. 
of , at point a’. 
of qa, at point d. 


nose skin. 

on upper skin of center cell. 
on rear cell skin. 

On bottom skin of center cell. 
on front web. 


on second web. 


Censider now the span wise forces acting on section pa of Figure 2, 


as shown in Figure 5. From equilibrium of spanwise forces: 


P — - P+ q a2 





-qdZ=0 
4 


dP 


az 








a 


ap 
where subscript ap indicates re- 
gion from point a to point b. 





-~M y 
x 
ol - al where a = Stress due to bemiling. 
x 
y = Distance from x axis. 
I. = Moment of inertia about x axis. 
M. = Bending moment about x axis. 
My ‘f we x 
- P= =. = SS 
Thus . dA r yaad Sc 
x x x 
ap ap 
where eo = moment of area from a to p about x axis. 
. dM 
But a V (i.e., shear force) 
x x 
VQ Q 
ence : eae (Be 
aZ i. x daz I 


Note: The second term in right side of above equation equals zero for 
non-tapered wing. 
Therefore: 


aad. + oe (2) 


Points a and a‘ as shown in Figure 2 are points adjacent to either side 
of webd-skin junction. Figure 4 shows the spanwise forces acting on section 





aa'p. 
eed? Jz From Figure 5: 
dz var 
| 
oP oo eg mies ct 2 8 
a a’ 1 aZ, I. 
va 
i Cae (2) 
x 








Considering spanwise forces acting on section a'p'' (Fig. 2 and ), ths 
equilibrium equation gives: 





ap 
P+ az : 
We BS haz] |, I 
a x 
or s 
VQ ty” 
qd, = 4, + 5 (3) 
x 


Figure 5 


Considering spanwise forces acting upon section p" g' (Fig. 2 and 5), 


the equilibrium equation gives: 





ap 
P+— dz 
az x 
ap VQiz. 3 
ag 45° (B Lee 
g OZ fF nie i 
Pg z 
from which 2 
VQ on, 
= + —2&- 
Tee 45 i 
x 
Figure ¢ 
x x x 
Ve rT 2 ae a 
and from equation 3, qa, 2a, tm + Peg + (B) 
g a I I ry I. 


Considering spanwise forces acting upon section gp"’ (Fig. 2 andé), the 


equilibrium equation gives: 





ap V2» 
Yr - = @ = = — 
Lo anu q, | azZ ap" : ae 


or from equation B, 


vq Va... 
ae 2 + eo (C) 
x x 





Figure 7 
Considering spanwise forces acting upon section g"p'Y (Fig. 2 andS ), 
the equilibrium equation gives: 


x 
VQ 
ep'y “x 


(D) 





Figure & 
Considering spanwise forces acting upon section p'Y d (Fig. 2 and 9 ), 


the equilibriwmma equation gives: 


x 
VQ ’ 
ap aP dp’ 
Pp > aZ = q ‘ BS = - 
az 7 4g * la eee i 
or 
x 
VQ 
a Vv 
4-7" (4) 
x 








From equations D and 4: 


x x x 
V V ¥ 
ang, RT ap'T 4. fed 
g” a I I. a I. 


And equation C becomes: 


x x 

vQ va. Vai, 
Wwe aso. ~~ eee: rt 

x x x 
or x x 

VQ VQ 
gp ' agfed 
= = + 

Wo * Ir 7 Ws % i i (9) 


Considering spanwise forces acting upon section dd'p"’ (Fig. 2 and 10), 


the equilibrium equation gives: 


aP 
- -- OZ, 
7 dz 





Figure 10 
From equations 5 and E: 
x x "3 
VQ it 9 VQ te | VQ 
a’ a I i I 
x 
or x 
¥Q 
agfedg 
= > = . r 
Tas ee I ) 


x 
Considering spanwise forces acting upon section d'p” (Fig. 2 andj1), 


the equilibrium equation gives: 





vQ. 
ie ene 
a'p’ x 


and from equation F 


x x 

VQ va. 
+ aes + er (6) 

x x 


c a! 





Figure 11 


Thus equations one through six represent the spanwise equilibrium equa- 
tions. These equations involve the unknowns 1? 1, » Wy? Wo? dp qd. 
qd? qs and qs° Four additional equations are now obtained by consider- 
ing the rate of twist of each cell and a moment equation. The moment 
equation will introduce a new unkmown (i.e., the shear center location). 





0 centroid 
8.c. = shear center 


Figure °° 


Referring to Figure )° and taking moments about point h: 


: = - dA _.- dA _ = 0 
(G) wer, 02 / a IB ot Sf; ¢.44 24 8 c2 | C3 
abc ag gted ac ed 


Substituting into equation G values of q> q,? 1.2 q2 Ve in terms of 


q? qo» and Qa" 





‘ 2 


x 
VQ tan Ve, iy 
aa -2{ la.,t SP laa -2 . = 
a 





‘ qa 
a I. fou i pied I. wei 





Veo 
a 
V + +2 mae 
{3 Xe.) ps a 
abc x 


x x x x 
vQ Won oe vQ va VQ, 
a! I I co a I ' a! I c3 
dc x x ga x x x 


x 





Dividing equation H by V and recalling qd? div and qs are constants: 











2 2 2 2 2 
(r+x = = at ~ ve’ fr ee + va’ A _- Sa A _- a Q™ dA 
8.c. Vou =v cl V ARL Vv ce Ve c5om ap n 
x jabec 
: 2q 
2 x 2 x 2 x a’ 
+7 teen 5 Qt foe Q e A 
a dg ce V Cc 
I ae Dec. I, afed dp qARy I ne agfedg 5 
2 x 2 x x 2 x 
eee Clee + —— Q dA = Q dA += ) dA 
2 d iL ORS 
I. ae ad'p" c I. Be gfed c3 te ea ag C3 = ee gp 5 
tA .+.A +A =A 
Bu a a °3 : (area of center cell) 
And Any - Aus = A. (area of rear cell) 
Hence, equation I becomes: 
re Xx ) = =_ ae, da - 2 @ A +A Sess A 
wc. o) =v os eoc.||«(O VAR I agfed | c2 c3] I, ‘gd 2 
2 x 2 x 2 x 2 x 
- = Q @a + — Ce an oe Q vist Tr QA 
Pcl of dp' 1 iI a. ce 
‘e abc i a “x ag = x | gfed a » ae es P 


The equations of the rate of twist for each cell give three additional 
equations. Note, the rate of twist for each cell is zero since the load 
is through the shear center. 


Front Cell: 
dQ ds ds 
2 = = —— = —— = {) 
04 || qt Wat 
abca adc ac 


19 


J 





1G i 


Center Cell: 








ao ds ds ds ds 
CGA carn = ee — . ams — = 0 
C Wa t tt Ve t tt 
acdg ac cd gd ag 
Rear Cell: 
ao | ds ds 
acs, | | a Q. oe cas a 
az t We t 
defg gfed gd 


Substituting into equations J, K, and L, the values of 


dp ? and q, in terms of q? qi. and q.. 


Ty? Ne? Wav? Ye? 


a 
Front Cell: 
ae 
V 
“ap ds Le ds 
ara arn -Q Ga ee 
a sE t a a I t 
abc x ac x 
or 
Taf as tet ( as af x as 2 x as. 
V t V . a Ape ii apt ete 
abca ac abc ac 


x x x 
VQ Q VQ ou, 
: ho a | re A 
ad I a! I. I t 
gd. ag 


al 





q q q Q 
ios). ae. 4 ds _agfedg | ds agfed | ds 
Vv t Vv t V t di t t 
ac acdga x ca x ed 
ut x Qe x ds 1 x d 
7 a Q tr7V¥ os Q rT) : = Q” ap = 
I ap’ t I d’p'* I gp’ t if i = 
ac ed ad ax 
Reer Cell 
= x x x 
VQ VQ vQ vQ 
a gf apt" 
see}. nay # yates TS a =o 
gfed >< ed. x x : x 
or 
a a Q™ 
t 
Vv. — Sige ds _ _agfed dg 1 = a aX as 
Vv t Vv t I t 1 ap'Yt 1 ag t 
efedg gd / ga gfed gd 





Thus equations I, II, III, and IV represent four equations with the 
four unknowns qo? qi qs and Ee Once these equations the desired 
shear flows can be obtained by use of equations 1 through 6. These 
equations then, apply to the solution of the shear flow problem involy- 
ing a structure of three cells, skin fully effective in bending, and the 
load acting through the shear center. 

To the shear flows previously calculated for the load through the 


shear center, those shear flows due to a torsional load must be added. 


Msc 
ie s.c. = Shear Center. 
Me é = Applied Torsional 
Figure 13 Moment about Shear 


Center. 
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From the previous spanwise equilibrium equations and the fact that 


d 
+ = 0 for pure torsion the following equations result: 


ied oa, 

* Wa i *o 4 ts 
ze a * a: 

10. q, = qs 

Be Wo * | 
i Wi 


Taking moments about shear center (Fig. 13) gives: 


(m) M +2. qnaA, + 2 Vito. > 2 
ac 


dA _-2 
9.C. al by c3 Lael 
abc fed 


ae deers 2 aA. = 0 
6 
cd oe gfed *R 


Substituting into equation ({m) the values of qe? Vey? Deo? ae and 
in terms of q 1 and qs gives: 


+ “pe -q- dA - 2 dA 
oF oc. : aa e ( qe a) ce Ts 
abc ac ag 


-2 (q .-q.)dA, - 2 q dA -2 aedée = 10 
: 4 ‘ ad 6 
gd <i cd i gfed 


But Ay . A, = AY (Area of nose cell) 


A + A +A +A = A Area of center cell 
c2 c3 ch c5 Cc ( ) 


A. -A, = A, (Area of rear cell) 
Hence: 


= A +2q A - 2a A 
Y) We ete Ma, Tan 
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Considering now the rate of twist in each cell and assuming that the 
rate of twist for each cell is the same gives: 
Nose Cell: 


ao ds ds 

2 =e = —_= a= => 

s) Se | Ay t qa t 
abca abc ac 





N 


Center Cell: 


(p) aca |S : es oo dg ds 
e c \az Wit 7 +t Wot 7 tt 
acdga ac da &. ag 
Rear Cell: 
add ds ds 
(aq) 2 rd = - a 7 
GAL \az Lae a tpt sf We +t 


Substituting into equations (m), (p), and (q), the values of Le? a,» q,9 


Wa? qo? and ah in terms of q qo and 14 gives: 


Nose Cell: 
de ds ds 
a2 GA, I ~ Qe t 7 Tye t 
abca abca ac 
Center Cell: 
de ( ds ds ds 


(VIT) Eley ‘| qos 1 qo af qa ae 
acdga acdga ac gd. 


Rear Cell: 


de ds 
a) 20k, td ~ ~ Ga oe Sas t 
gdefe2 efedg ed 


Thus, equations (V), (VI), (VII), and (VIII), represent four equations 


de 
and four unknowns, 1? qo q., and rd - Once these equations are 


a’ 


sOlved the shear flows due to a torsional load about the shear center 
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may be determined by use of equations (7) through (12). 
Summarizing the equations for the determination of the shear 
flows in a three cell beam with skin fully effective in bending gives: 


A. Equations for Vertical Load Through Shear Center: 


2q A 2q A 2q A 





a N a'c akR.e2e x 
I + = - te st A +A 
> ie.c.! v v VT “agfea | o2 °3 | 
a x 2 x x 2 De 
+-_- .- = = _- = 
r eta E cor Ma’ " Gee cl IT Sap AR 
abc ag gfed 
2 x 
* 7 Sarpy Aco 
dc 
(II) * de | a! de i o> ds Jl oX a8 9 
t et 
i abca x ac ax abc *P t ac *P 
(III) wa ae Ua! de jaf ae x ds 
Vv t V t V t agfedg t 
ac acdga ed cd 
ee 
agfed ds ae ot de lt aX de 
a t I ap' t I d'p'¥ + 
: ed ac cd 
1 x de 1 x ds 
_~ —_-=— = 0 
I Sept t I Sa tp" t 
ed. ag 
a 
(Iv) Sa as _“e'( as “agfea ( ae 1 x as 
v te * t I t ap'Y t. 
gfedg gd ed efed 
Lie & | « 
I. ag t 
gd 





x 
2) a A+ ‘a. 


(2) oy ae cel. e 
(3) ae, * 


ote yy 


“ fed va wae 
, —2gfed | _“gp'" 
I = 


(>) ae = Tees 


vQ" vQ. 
Cc a’ I I 


B. Equations for Torsional Load: 
V =2q A +2qg A - 20 A 
(V) ite. Qo c qs R q, N 


de 
(VI) 2ca.|— 








R 


(VII) 26a 





(VIII) 26a, 





(7) Q,= 4, 


(8) Wy “a i qe 
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(11) MP, 4 > Oy 
(12) Cea a 


do 


az, 


do 


(13) 








abca a acdga gdefg 





Referring to the above summary of equations, the resulting shear 
flows for pure torsion are of constant value between any two successive 
etructural junctions, while the shear flows due te bending vary between 
such junction pointe. The maximum and minimum values of these shear 
flews can of course be obtained by differentiating the shear flow ex- 
pression with respect to the parameter defining the amount of area of 
the appropriate section. 

The procedure employed in the analysis of the subject three cell 
beam is applicable to a 'm' cell beam (fer n > 2) where, 


3n - 2 = number of unknowns (1.¢., shear flows and shear center 


location). 


n = number of equations available from the rate of twist of 
each cell. 
1 = moment equation. 


Thus: 


(3n - 2) - (mn - 1) = 2n - 3 = number of points at which the shear 
value is assumed. 


Previously it was shown that when considering the principal axes 
x 
V 
She - 
az} I 
ap x 


then 





If the axes are not principe] axes and in addition the load is acting 


in any given direction then 
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ap y 
MI -~M I M I ~ “I 
= : xdA + iI : yaa 
y ap J xy ap 
- MI vi I - MI ‘ 
Ef me + —_27___,.*_F 
- a ro c= ee 
x ty xy xy 
1 x x J 
Pp = M (I - I +M (I =F 
II - ite | “i att bp) uA “cers Ay| 
x¥yr | Cexcy 
ap 1 y Ae 
a = - + V - I 
(IX) az Il -I 2,58, it = >) ie ee) 
x y xy 
Te ge 1G. T spe - ee 
‘aK a xy ap v8?) oy a| xy a x a 
x daZ II - y az II - 
xJ xy xy xy 


Thus the algebra now involved is very cumbersum, and a practical 
analysis would require some assumptions or approximations wherein 
this difficulty could be overcome. For example, in the procedure 
employed herein it was assumed the wing or beam had no taper, hence, 
the lest two terms of equation (IX) each equals zero. Also principal 
axes were used, thus te equals zero, and only a vertical load (1.e., 
perpendicular to the x axis) was considered. These qualifications 


reduce equation (IX) to 


vQ 
ap _ 8p 
af oS 
ap x 
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Alse worthy of note, is that the shear center location as derived 
herein is different than that location as obtained from the analysis 
Where the skin is considered either nen-effective or pertially effec- 
tive in bending. 
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PART II 


STRESS ANALYSIS AND STATIC TEST OF A MODEL WING 
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PART II 


STRESS ANALYSIS AKRD STATIC TEST OF A MODEL WING 


This part of the report includes the stress analysis and static 
| test ef a three cell model wing. The wing was constructed from 

.O72-24ST aluminum sheet, with six foot span, twelve inch chord, and 
circular arc airfoil shape. As far as the stress analysis is con- 
cerned, those equations developed in Part I could be applied directly, 
however, due to the symmetry of the wing, the shear flew equations 
are somewhat simpler than those of Fart I, hence, a separate stress 
analysis for the wing using the previous procedure is included. The 
initial sections of Part II are deroted to this stress analysis and 
the actual static test follows. 


Geometric Properties of Model Wing 


Due to symmetry about chord line 
only upper half of wing need be 
considered. 


(y+a)* + x” = Re 


atx = 0,y = .64 t= .072" 
R = 28.445" 
6h+e = R c= 12" 
a=R- .64 = 27.805" a = 27.805” 
sin y= fe ee R* = 809.118 sq.in. 
RO 26.445 a? = 773.118 ag. in. 


= .9175 7 = 1.3585 rad. 


sin y = .9171) 
y = 1.3585 rad. aR = 790.9122 








Moment of Inertia (I): 


We 
aul e oan 
ar egicin | i = ta 
7 


Wo 

2 _ 2\2 

ea | (R sin g - a)*~ Rta 9 
1; 
2 2 
2 unt | EB - + cin of + 248 cos J + 079 
1. 3585 
= 8.192 {404.559 [1.570796 - 1.3585] - 202.279 [0 - sin 2.7170] 
+ 790.9132x 2 [0 - cos 1.3565] + 773.118 x .212296 | 
= 8.192 {404.554 x .212296 + 202.279 x .411949 + 158.18264 x 
[-.21070| + 773.118 x .212296} 

= 8.192 } 85.885 + 83. 32863 - 333.29082 + 164.1298} 

= 8.192 x .05261 = .4309895 in.* 

_ 4 
ee - 4309895 in. 

- bh> 072 K .96° © ; 

ie = 2 ae zZ = .012 x 875136 = .0105016 
I. Soe: raat] = 4 aE + .6X O72 x .475 

= .0097648 ) 

? 4 
Ie stay = MLMg An. 
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Geometric Properties of Model Wing: 


a) Note: The rivet line for the spar 
ie a is 4 inches from trailing edge. The 
‘= Web location is 3.58 inches from 
< trailing edge. For purposee of analy- 
a ee sis the spar loading is assumed to 
(ae ail transfer to the skin at some point 
a between the rivet line and web line 
(e.g., 3.79" from t.e.). 
Sin! 2 o> Senora aie. ° . ae 211 
~ RO 28.NN5 “ R 28.445 
W = .0778 rad. 8 = .2126 rad. 
Lt = .64 


SO = RcosW = 28.445 .99698 = 28.3591 
av = st = .5541”" 


y 
uw | RdOQ = R(Y - 0) 
La! 
— [= 
Y 
“a = 7.604" 





° 
.@) 
@ 

<= 
i 


99698 cos@ = .9775 


28.445  .0778 = 2.215" 


28.445 (.2116 - .0778) = 3.802" 


1 
Fad 
Ww 

1 
a 
il 


= 4.430" 
La’ 
ce 
=O | 
= A : a = a or Area 'sOL' = OTB * 809. S 31.47459 
Cc R 


2) 





Area "aos" = 4X 28.3591 X 2.21 = 31.3368 


2.21 *% .55 + 31.4759 - 31.3368 = 1.3522 sq.in. 


rl~- 
> 
tt 


A = 5.4128 sq. in. 


4 A = Area "boL” - 2 A, > "bet" 


n 
Area _'boL' 6 area "beL" = 22226.% 809-118 _ g¢. 00924 
R* eTT 2 
Area "bot" = $X 6 Xa = 3 X 27.805 = 83.415 


BA = 86.0094 - 1.3532 - 83.418 = 1.2h07 
A= 2.4814 sq. in. 


Calculation of Moment of Areas: (Refer to Fig. 6) 


9 9 
QF oy a8: Qu = y Rt dO = Rt (R cos G@ - a) dO 
tet P ) O 


Rt [R sin o - ad} 


YW =.0778 
| a ds -| R*t [R sin 0 - ad| do 
La' 0 


; .0778 
R°t . R cos 9 - 4 
‘ 0 


| 


809.118 x .072 { 28. uu5( 99698 ay 
_ 27.805 ~— | 


2 


fl 


58.2569 { 0859039 - 0841465} 


- L0237995 


2k 





9 70 
x x 
| Geds: “9 | Rt ydO = Rt (R cos 9 - a) dO 
ap ap 
ab y y 


Rt [R(ein 0 - sinW) + aly - 9) | 


p 
| Q as = Rae [R(ain 0 - sinW) +a (W- 0) | de 
ab 





ap 
iy all 
- R°t [- R cos 0 - RO sin +a 0 - ee 
y 
- R*t | - R(cos @ - cos) - R (8 -Y) sin’ +awP 
2 2 


58 .25649 {- R(-.01948) - R(.1348)x.0777 
2 2 
+a x .0778 x .21z6 - S1-9178)" _ 9( 2126)" 


58.25649 {.5541086 - .29793179 + .459902K8 
- .0B41496 - .6283531} 


98.25649 X .00357659 


= 208359579 


h=.5541 2 2 
t t 
j Q7 de: Qe yt ay = HY. ES 5 (n? - y*) 
&C y 


h 
ie , ads = 2 | 
ac P 0 
= ef ae 


a 
7 2 x = x .072 = ,00816592 


ML h 
F (h* - y*) ay = t sy - rr 
0 
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3 3 


2) 











Fig. 


>= 
oe) 
me 
tt 
Qu 
@ 
1] 


a‘L 


W=.0778 y 

xX 

Q | | oat | Rt(R cos 0 - a) ad 
9) O 


.0778 
Rt E sin 6 - wo | 
QO 


2.0480K [28.445 x .O777 - 27.805 x .0778 | 


1 


2.04804 [ 2.2101765 - 2.163229] 


. 096147456 


09614745 x 1.10 


. 10576219 


Geometric Properties of Wing 
(Strain Gage Locations) 


- on 
968 
14.9 (138-5 


425 





x 
ae 
6.03! ws gt it _|| 6 
344 | y | | 


eS 
2 Planform of wing showing strain gauge locations: 
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Load Point 








Referring to Fig. 3, the strain gauge locations are those ac- 
tually used in the static test. They are introduced here in order 
to analytically determine the shear stresses at these locations. 


Stress Analysis for Vertical Load Through Shear Center 
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eV x V x 
4 ds + = = = 
| + ae ds TL Se ds <| (Qo + oe ds = 0 
ab ac jab ac 


2q_(3.802) + q (1-10) + = ( 2083596) : = ( 0081659 + ,10576219) = 0 


x x 


8.704 q += [ .m67ige - .1239281] = 0 
x 


8.70 q. + .3027911 = = 0 
x 
v 
fe - O01) 
x 
a = -.0788 V 
x 
VQ 
a 





Y 

'L 

y on | Rt (R cos 9 - a) do = Rt[R siny - ay | 
O 


2.04804 [28.445 x .0777 - 27.805 * .0778] 


2.04804 [ 2.210176 - 2.163229 | 


= .0469475 
O46 
1a "Sag 
q . = .1065 ¥ 
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Strain Gage C. 
aig =e Se. ganee 


R 28.445 — 
fk 3 75> . a= O48 rad. 
| iy * Tigre - on eee 






an = Rt[R sin o - ad] 


QF = 2.08604 [28.445 x .0h83 - 27.805 x .OK8K | 
= 2.0h804 [1.37389 - 1.343762 | 
Figure 5 
= 2.04804 *< .03013 
; vor = .06171 
| Rie 
gage c x 
Te = 13977 V 
Eo , Fe 1-2 


p b.75 
ak aS ae 
x 
Be, =. 
yor ~ BB NAS -1669 
“oo oo 


Rt [R cos 0 - a 40 


= > 


+ Rt | R(sin e - sin) - apt aly] 

2.04804 [28.445(.1679 - -O777) 
- a(.1679 - .0778)| 

2.0480k [2.565739 - 2.5052305 | 

2.04804 x .06051 = .1239269 





eg 





Hence, the shear stresses as analytically determined at strain 


gauge locations B and C, due to a vertical load through the shear 
center of the wing are: 


At lecation B: a 2.804 VY p.s.i. 


At location C: Te 1.941 V p.s.i. 


it 


where V = applied load. 


Stress Analysis for Moment about Shear Center 





Due to symmetry q = - a, (d) 


Figure 7 
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Taking Moments about Point 0: 
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z dA +2 dA +2 dA _ -2 dA =" 
"B= cl eee tl Ay op te cs mee 
La’ abc ac Je'e 
q dA + da. + ( ) You 1. 22 
a’ cl aq nl Qos a. c2 Qo c% Ok 
La' abc ac a'k 
dA + + aa _| + aa ae | see 
Gas cl ce C3 q, nl ce} 4& 
La ' ac c'k abc ac 


Jk 
+q A == = 
>) qa. A, +4, 4,7], where A. 4 area of center cell 


area of nose cell. 
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Twist of Nose Cell: 





Figure & 
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aca [42 2 co 4 t 
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az 2GA t q t ae qe 


abc ac ac 








Twist of Center Cell: 











2GA Se Ss + Es + Ls + S 
c| az tt Wt tt Wt 
foe} ac cd dg 
2 2 
= 7 ds + — 
| Sepas | ay, “5 
ga ac 
2 2 
7 an Tee a0 (a, - a) ds 
FP, a 
= 42 “Ue 7 , ae | ds 
ga AC 
ao 
— = + - 
ga ac ac 
| ae " 
ee nose a center 
cell cell 


1 ut 
an | 4, (7-604 + 1.10) - a.% 1.10 | = rv | a, «(4-430 + 1,10) ee ae 1.10] 


4128 


ea 511 (8.704 co 1.10 a.) = 5.530 a 1.10 q. 


W.THESS a - - 59X65 a., = 9-520 4., - 1-104. 
5.84655 qs. + 6.12986 qo 
(6) 2» 1.0484 qi 


From eqt. (5) 22828 OF is 2.4814 x 1.O48kq - * 


5- 30789 4, = 


. Fly 


q., = -O471 T Teh, One® 
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q. = 1.0484 q , = 1.0484 X .0471 


q_ = .04938 T T, = -6858 1 
» 
Hence, the shear stresses as analytically determined at strain 
gauge locations B and C, due to a moment about the shear center of 
the wing are: 


At location B: Ve .6858 T p.s.i. 


.6542 T p.s.i. 


At location C: ‘ie 


where T = applied moment. 
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EQUIPMENT AKD PROCEDURE 


The geometry and material of the model wing has been discussed 
previously. An attempt was made to test the wing in pure bending 
(i.e., no tersion) and in pure torsion. Six SR-4 strain gages were 
mounted on the wing and an additional one installed upon a small sheet 
of 2kST aluminum for temperature correction purposes. (See Fig. 3 for 
strain gage locations.) The rig by which the bending test was con- 
ducted is shown in Fig. 9. Referring to Fig. 9, the load was applied 
by adjusting the turnbuckle which appears on the right. The amount 
of load was determined by a previously calibrated dynamometer installed 
just below the turnbuckle. Strain gage readings were obtained at 
several load values, varying from zero to 437 pounds. (See Tables I 
and II at end of this report.) 

Following the bending test the wing was tested in torsion. The 
rig used here is shown in Figs. 10 and 1l. Referring to Fig. 10, 
the load was applied by a hydraulic jack and amount of load deter- 
mined by use of the dynamometer used in the bending test. Not ap- 
parent from Figs. 10 and 11 1s the balance scheme used, wherein a 
cable was attached to the free end of the wing and loaded with sand 
bags such that the wing was not subjected to bending due to its own 
Weight. Strain gage readings were obtained at several torsional 
load values, varying from zero to over 6000 in.lbs. (See Table III 
at the end of this report.) 
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‘ig-10. Torsion Test 
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RESULTS AND DISCUSSION 

rom the resulting test data (i.e , Tables I, II, and DI) Mohr 
stress circles were drawn for each strain gage, corresponding to 
V = 400 pounds in the bending test and T = 6230 inch pounds in the 
torsion test. This was done somewhat as a preliminary check (1.e., 
a comparison between experimental shear stress values with those 
calculated previously). Considerable disagreement in results was 
indicated at thie point. It appeared that perhaps the bending test 
actually introduced some torsional moment upon the wing and similarly 
the torsion test introduced some bending upon the wing. Having suf- 
ficient data, it was possible to attempt to separate the torsional 
shear stresses from the bending shear stresses in both the bending 
and torsion tests. This was done and the experimental results were 
compensated for accordingly. However, there was still considerable 
disagreement between the experimental and analytical results. The 
stress circles are not included in this report since their only ap- 
parent value was to show considerable disagreement between the ex- 
perimentel and analytical results. 

In view of this disagreement in results, the next logical ap- 
proach appeared to be an attempt at determining the reliability of 
the experimental strains. Figure 12 is a plot of load vs. strain 
from strain gages C, B, and E. From this plot it is apparent that 
there was no over-all linearity of stress-load variation as pre- 
dicted by Hookes Law. The causes leading to this non-linearity be- 
tween stress and load might be attributed to the inherent random 
error in the strain determination, instability of the structure 
wherein the slippage at the riveted joints causes portions of the 


structure to load or unload in some unpredictable fashion, error 
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In discussion each of the above sources of error, the results 
of Figure 12 indicates inherent random error of the strain 
determinations was of minor importance, since the irregularities 
in the strain are of much larger magnitude than those attributable 
to random error, where it is known that the strain-meter readings 
as such, are at least accurate to within plus or minus one-half 


microinch. Error in the load determination may be ruled out since 





Such errors would need be in order of magnitude of plus fifty pounds 
at some points and minus fifty pounds at other points. 

Errors due to scale factor of the model and limitations of strain 
Gages appear as possibly large, when considering that the magnitude 
of shear stresses is quite low and the shear stress values so critic- 
ally dependent upon the precise location of the strain gages. For 
example, in bending with the applied load equal to approximately 
four hundred pounds the shear stress in the skin over the center 
cell of the wing varies from zero at the center of the panel to ap- 
proximately eight hundred pounds per square inch at the web skin 
junction. Thus the shear stress varies from zero to eight hundred 
pounds per square inch over a distance of about two inches, hence, 
the requirement for precise location of the strain gages. In this 
respect it should be noted that when using a strain rosette the 
strain readings are assumed to represent those strains at the inter- 
section of the gage axes on the strain gage (Fig. 13). 

Figure 13 is a full scale diagrammatic representation of the 
strain gages used in this experiment. Consider any gage axis (e.g., 
gage axis 3) and suppose that when mounted on a structure the strains 


in the 3 direction have considerabls variation over the gage length. 





, 
Fig. 13 

The strain gage reading in the 3 direction would indicate, approxi- 
mately, an average strain over the gage length, and as previously 
stated, it is assumed that the resulting strain is representative 
of the strain at the gage axes intersection However, it is quite 
possible that a truer location of the resulting strain would be 
somewhat to the right of the gage axes intersection point, especially 
if the lead in wires to the gage have a portion glued to the sur- 
face of the structure. From close observance of Figure 13, it ap- 
pears that the true location of the resulting strain may be as much 
as @ querter of an inch to the right of the gage axis intersection 
point. In usual applications of strain gages this error would be 
negligible, however, as previously mention, in the subject test the 
strains might vary considerably over the distance of onle one-quarter 
inch. This error, if present, cannot be compensated for, since if 
the true strain locations along the axes are not at the intersection 
of the axes, then data taken from the gage readings would be used for 
constructing one strain circle to represent three different points 
on the structure, which is meaningless. But such error can be dis- 
counted as the main source of error in this experiment, since such 
an error would not cause the abrupt changes in slope of the strain 
vs. load curves, in fact, the *rror would be linear with load. Also, 
there is present the fact that such error could not account for the 
discrepancies in the torsion test, since here the shear flows are 
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constant over any integral section of a cell (by integral section 
is meant the skin as opposed to the webs). 

The error due to the local effect of the rivets could certainly 
be prominent for some of the gages but would not explain the discrep- 
ancies shown by gage B which was located at the middle of the center 
cell skin and at maximum distance from a rivet line. 

The instability of the structure, wherein slippage at the riveted 
Joints causes portions of the structure to load and unload in some 
unpredictable fashion remains as the possibility of the cause for the 
large variations in stress-load slope. A further study of Figure le 
makes this appear as a reasonable conclusion. For example, note the 
strains at the loads of 425, 363, and 250 pounds. Ideally, these 
strain values would form a straight line for each gage axis, but if 
a straight line is drawn between the 250 and 425 pounds strain values 
the 263 pounds strain value does not lay upon this straight line for 
any of the gage axes. The 363 pounds strain values do though, show 
a uniform trend upon considering that gages G, B, and C have the 
same gage axis direction orientation on the wing, and the 363 pounds 
strain values all lay to the left of the aforementioned straight line 
for gage axes B,, C,, G, and Bo, Co, Gp, and to the right for gage 
axes Bs, Cs, and Gs. Here again is evidence that the error is not 
caused by either random error or error in say any one set of gage 
readings, but rather that the structure was behaving as the data 
indicates, such behavior being induced by slippage of the riveted 
joints. 

The slippage of the riveted joints could be expected to not be ~- 
continuous with load, but intermittent, such that portions of the 
structure would "lock" over certain load ranges, slip and lock 
again. Thus, it is probable that over some loading range, during 
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which the structure remained in a locked condition, there would be 
the expected stress-load linearity, and such is the case for the load 
range of 128 to 250 pounds. Hence, it appeared that for some load 
range the test data might be representative of the structure without 
slippage of joints. To determine this range the plot of Figure 14 
was constructed. The analytical values of strain as shown there are 
included in Table IV. From Figure 14, it appeared that the mst 
probable load range in which the experimental strains would repre- 
sentative of the actual wing without slippage was between 166 and 250 
pounds. Even in this range the experimental and analytical stress- 
load slopes are not in agreement but this can be explained from the 
fact that even though the test was intended for pure bending there 
was somes torsion entering, as indicated by the strains from gages 
Gand F. The next step was to determine if the amount of torsion 
entering would compensate for the difference between experimental 
and analytical stress-load slopes. This was done from analysis of 
Figure 15. 

From Figure 15 note the experimental load-strain line for gage C. 
This was corrected for torsion by subtracting those strains due to 
torsion as indicated by the shearing strains from gage G. The re- 
sulting load-strain line showed close agreement with the analytically 
derived line. However, the same corrections when applied to the 
strains of gages B and D do not give agreement between the experi- 
mental and analytical results, hence, the indication that slippage 


occurred over this load range also. 
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LOAD vs SHEAR STRAIN 
EXPERIMENTAL & ANALYTICAL RESULTS 
Figure 15 





¥¥e 2.0 a. a 56 
SHEAR STRAIN (x10° ) 


© 4.0 "7 








CONCLUSIONS AND RECOMMENDATIONS 


The results of the static test of the wing are of little value 





but do indicate some recommendations for such a test. For such a 
test the scale of the model should be larger in order that the "life 
size" rivets do not predominate upon the stress distribution of the 
model, and also to reduce the preciseness requirement of the strain 
locations, as previously discussed. Also care should be taken to 
provide rigid joints, so that the slippage effects as evidenced by 
this experiment would be minimized. In fact, it appears that a dif- 
ferent method of joining the skins at the leading and trailing edges 
would be preferred to riveting. 

Another recommendation is to have the spar flanges as small as 
minimum rivet edge distance would permit. This conclusion arises 
a condition whereby the load transfer between spar and wing skin 
does not necessarily take place at the rivet line. For Saiie. 
suppose that the skin-spar force relation is such that the two sur- 
faces are being pressed upon each other, as opposed to being drawn 
apart. The load transfer between the two members would then be 
somewhat removed from the rivet line and distributed near the web 
end of the flange, hence, the desirability of keeping the flange 
width small. 

The above recommendations are in no way intended as recommenda- 
tions for wing construction as such, but are suggestions whereby 
closer agreement between analytical and experimental resulte might 


be realized. 


AG, 








Strain Gage 
Side in 


Strain Gage Side in 


Strain Gage Side 


in Compression 


Strain Gage Side in 


Tension 
ey 


Tension 
= 


Compression 


= Dynamometer Strain Readings. 


RABD 
Mm OKO AN @ 


Ww 
ON 
Wi so Wn 


> 
FO 
oo 


g 


ae 
as 
OO 0 F 


: R 


RSPB 
oO OT © 


BBA 
FOoWr 
Oon oF 


OCOOF 


SIAS 
ore? ta 


689.7 


619.6 


695.2 
697.9 
704.9 
707.3 
711.2 
720 5 
T27.5 


691.8 
678.8 
671.5 
656.9 


485.0 
48h 8 
488 . 3 
490.2 
495.9 
496.6 
499.0 
505.7 
510.3 


487.8 
477 Ju 
472.3 
461.9 


z 
2 


508.8 
510.3 
514.0 
516.6 
520.2 
9ee.9 
529.8 
Boi «2 
Suy. 3 


905.2 
4O7 .6 
490.8 
477.0 


472.7 
472.1 
471.0 
470.0 
472.4 
470.7 
466.1 
465.8 
463.6 


479.8 
476.9 
480.2 
LBL. 3 


3 


463.0 
461.9 
459.0 
457.1 
455.0 
AS 1 
450.5 
ALT 4 
LhS, 3 


466.5 
471 5 
479.9 
Ye. 1 


439.9 
ban. 1 
455-5 
460.1 
468.4 
475.5 
4OL.1 
D077. 0 
221.7 


440. 3 
412.1 
391 3 
DILae 


619.3 
625.2 
636.6 


658. 3 
665.0 
683.3 
705.1 
723.8 


ag 
NM OO 


” 


2 
2 


ON 
OV © 


al 


Strain Readings in Microinches 


BENDING TEST 


B, C, E, F, and G = Strain Gage Readings. 
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STRAIN GAGES (all atrain readings 








Dyn. ibs. B 

1 a 2 2 2 3 
362.2 O 
362. 12.0 0 698.7 0 485.0 0 472.7 O 
364.2 41.0 29.0 691.6 1.9 484.8 -.2 472.1 - .6 
365.6 69.0 57.0 695.2 5.5 Youn 2.5 te eee 
366.9 95.0 83.0 697.9 8.2 490.2 5.2 470.0 -2.7 
369.1 140.0 128.0 704.9 15.2 495.9 10.9 472.1 - .6 
372.0 178.0 166.0 707.3 17.6 496.6 11.6 470.7 -2.0 
375.2 262.0 250.0 711.2 21.5 499.0 14.0 466.1 -6.6 
380.8 375.0 363.0 720.5 30.8 505.7 20.7 465.8 -6.9 
364.0 437.0 425.0 727.5 37.8 510.3 25.3 463.6 -9.1 
Dyn. lbs. C 

1 3 2 2 B, 3 
362.2 ©) 
362.8 12.0 0 439.9 0 417.7 0) 394.8 @) 
364.2 412.0 29.0 4hb.2 4.2 419.5 1.8 392.8 -2.0 
365.6 69.0 57.0 453.5 13.6 423.7 6.0 389.2 -5.6 
366.9 95.0 83.0 460.1 20.2 426.8 9.1 387.3 -7.5 
369.1 140.0 128.0 468.4 28.5 431.0 13.3 384.9 -9.9 
371.0 178.0 166.0 475.5 35.6 434.9 17.2 384.5 -10.3 
Bfo-e 262.0 250.0 Hola. 51.2 WhO.1 o8.h 377.9 -16.9 
380.8 375.0 363.0 507.0 67.1 449.3 31.6 371.6 -23.2 
784.0 437.0 425.0 522.7 81.8 456.8 39.1 366.7 -28.1 
Dyn. lbs. D 

1 i 2 2 2 3 
362.2 0) 
362.8 12.0 0 734.9 0 429.9 0 2h2.0 0 
364.2 41.0 29.0 734.4 -.5 430.4 21.4 243.8 sé 
Bop-0 69,0 57.0 733.5 -1.6 433.0 4.Q 2h7.6 5.6 
BO6.9 995705 65.0 731.6 -3.5 434.5 5.5 2h9.9 7.9 
369.1 140.0 128.0 730.3 -4.6 436.6 7.6 253.2 11.2 
371.0 178-0 166.0 729.6 -5.3 438.8 9.8 255.8 13.8 
375.2 262.0 250.0 726.7 -8.2 443.2 14.2 262.4 20.4 
380.8 375.0 363.0 724.9 -10.0 448.3 19.3 269.3 27.3 
384.0 437.0 425.0 722.3 -12.6 452.4 23.4 276.2 34.2 
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P = Dynamometer Strain Readings 






B,C, D, D, F, and G = Strain Gage Readings 









P B 


1 2 ° 1 2 


354.5 691.3 494.7 480.7 435.6 416.8 399.8 732.9 423.7 237.7 
355.9 690.1 493.7 480.0 435.0 416.7 399.2 733.2 423.7 237.7 
357.3 690.2 499.5 480.0 435.4 hO7.6 397.5 733.8 417.7 236.8 
361.7 690.2 512.1 479.5 436.9 8521.2 391.3 735.3 405.0 235.3 
365.6 690.4 524.3 479.0 437.9 474.1 387.6 736.8 391.8 234.7 
368.3 692.3 532.1 479.3 438.8 487.3 384.1 737.7 385.0 234.2 






C 






1 

354.5 680.7 508.12 465.6 615.5 358.9 337.7 648.9 470.9 450.6 
355.9 680.5 508.7 466.1 613.9 358.8 337.9 648.9 470.3 450.4 
357-3 680.9 514.4 464.7 624.3 371.3 337.5 647.7 482.3 450.8 
361.7 681.7 529.0 463.5 614.5 396.5 333.8 646.7 516.3 449.9 
365.6 682.8 544.8 463.6 613.9 417.2 307.4 645.8 528.1 446.6 
368.3 683.9 554.1 464.4 634.7 4207.8 3201.6 645.7 539.9 443.5 


Strain Readings in Microinches Table III 
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Table IV 


Shearing Strains for Bending Test 
V = Applied load in pounds. 





. Gage C 
V 2 2 10° -( y+ 3) 10? Pe 10° Peae 10° 
0 ©) 0 0 @) 
29 3.6 2.2 Teast 1.41 
57 12.0 8.0 4.0 2.90 
oo 18.2 12.7 5.5 4.20 
128 26.6 18.6 8.0 6.60 
166 34.4 25.3 Gul 8.50 
250 48.8 3h. 3 14.5 12.81 
363 73.2 43.9 29.3 SAG. 
425 78.2 Sp fi eh.5 21.75 
Gage B 
V2. 10° =( ae) 10° —_ 10° le 10° 
) 0) 0 0 0 

29 -.4 1.3 = ihe 76 2.10 
57 6.6 3.8 2.8 4.20 
83 10.4 5.5 4.9 6.22 
128 21.8 14.6 7.2 9.55 
166 23.2 15.6 7.6 12.25 
250 28.0 W.9 ilo sat 18.50 
363 41.4 23.9 LS 26.90 
425 50.6 28.7 21.9 e250 
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Stress analysis and static test of 
a three cell beam. 
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